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ABSTRACT 

Using third-order perturbation theory, we derive a relation between the mean diver- 
gence of the peculiar velocity given density and the density itself. Our calculations 
assume Gaussian initial conditions and are valid for Gaussian filtering of the evolved 
density and velocity fields. The mean velocity divergence turns out to be a third-order 
polynomial in the density contrast. We test the power spectrum dependence of the co- 
efficients of the polynomial for scale- free and standard CDM spectra and find it rather 
weak. Over scales larger than about 5/i _1 Mpc, the scatter in the relation is small 
compared to that introduced by random errors in the observed density and velocity 
fields. The relation can be useful for recovering the peculiar velocity from the associ- 
ated density field, and also for non-linear analyses of the anisotropies of structure in 
redshift surveys. 

Key words: cosmology: theory - galaxies: clustering - galaxies: formation - large- 
scale structure of the Universe 



1 INTRODUCTION 



In the gravitational instability paradigm for the formation 
of structure in the Universe, the peculiar motions (devia- 
tions from the Hubble flow) of galaxies are tightly related to 
. ^ the mass distribution in the Universe. Quantitative relations 
, between the peculiar velocity field, v, and the mass density 
5_l ■ contrast field, 5 — p/pb — 1, where pb is the background 
density, can be obtained from the equations of motion of a 
collisionless self-gravitating system. In the linear regime, the 
relation between the density and the velocity fields is 



s(x) = -f(ny 1 v -v( x ). 



(i) 



where /(fi) ~ Q 6 (e.g. Peebles 1980) and we express dis- 
tances in units of km s - . Given an assumed value for Q, the 
above relation can be used to reconstruct the mass density 
field from the large-scale velocity field. The comparison of 
the reconstructed mass field with the observed large-scale 
galaxy density field serves as a method for estimating fl 
(Dekel et al. 1993) and as a test of the gravitational insta- 
bility hypothesis (but see Babul et al. 1994). The linear the- 
ory relation is applicable only when the density fluctuations 
are small compared to unity. However, the observed den- 
sity fluctuations from current redshift surveys (e.g. Fisher 
et al. 1994) and from the potent reconstruction of density 
fields, slightly exceed the regime of applicability of linear 
theory. For example, the density contrast in regions like the 
Great attractor is about unity even when smoothed over 



scales of the order of 1000 km s" 1 (Dekel et al. 1993). Fu- 
ture redshift surveys and catalogs of peculiar velocities will 
provide reliable estimates of density and velocity fields on 
scales where nonlinear effects are certainly not negligible 
and need to be incorporated in analyzing the data. Various 
nonlinear relations between the velocity and the associated 
density field have been developed. One approach is to as- 
sume phenomenological parametric forms of these relations, 
which are to be calibrated with N-body simulations (Dekel 
1994, Ganon et al. 1998). Another complementary approach 
is to analytically derive these relations based on various ap- 
proximations to nonlinear dynamics. For the purpose of de- 
riving the density from the observed velocities, the potent 
algorithm for example uses a nonlinear approximation based 
on the Zel'dovich approximation (Nusser et al. 1991). The 
inverse problem, in which one wishes to recover the velocity 
from a given density field is yet unsolved in the Zel'dovich 
approximation as it involves a set of nonlinear differential 
equations, which do not have analytic solutions. This paper 
aims at deriving a nonlinear relation which can be used to re- 
cover the velocity field from the density field. This is relevant 
to comparisons of observed velocities of galaxies, using for 
example the Tully-Fisher relation, with predicted velocity 
fields from the density field estimated from redshift surveys 
(Strauss & Davis 1988, Yahil 1988, Kaiser et al. 1991, Hud- 
son 1994, Davis, Nusser & Willick 1996, Willick & Strauss 
1998). Moreover, such a relation can be useful for estimating 
SI from the distortions of clustering in redshift space (Kaiser 
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1987, Fisher & Nusser 1996, Taylor & Hamilton 1996, Hamil- 
ton 1997) when measured, for example, from the anisotropies 
of the correlation functions or, equivalently, the power spec- 
tra. This is particularly important for measurements of red- 
shift distortions in the Sloan Digital Sky Survey (Gunn & 
Knapp 1993) and the Anglo- Australian 2dF galaxy survey. 
Nonlinear analysis of these data might also aim at breaking 
the degeneracy between £1 and bias (Chodorowski & Lokas 
1997; hereafter Paper I, Bernardeau, Chodorowski & Lokas, 
1998). 

One way to go in order to derive velocity from density 
is to assume that the divergence of the velocity field at any 
point in space can be approximated by an expansion in terms 
of the density contrast at that point. The coefficients of this 
expansion can then be determined using either N-body simu- 
lations (Nusser et al. 1991, Mancinelli et al. 1994) or pertur- 
bation theory (Bernardeau 1992) under the assumption of 
Gaussian initial conditions. The advantage of this approach 
is that it provides the velocity divergence directly in terms 
of the density. For irrotational flows, the velocity field can 
readily be recovered from its divergence given some bound- 
ary conditions at large distances. Indeed, if we define?] 



9(x) = -/(n) _1 V • v(x), 
then the velocity field is simply 

m) 



v(x) 



4tt 



d 3 x'6(x')- 



(2) 



(3) 



One caveat to this approach is that the velocity divergence is 
not determined uniquely by the density field (Chodorowski 
1997, Mancinelli & Yahil 1995, Catelan et al. 1995) and the 
scatter around the derived relation is expected to propagate 
into errors in the velocity field. However, in the weakly non- 
linear regime, the values of 5(x) and 9(x) are still strongly 
correlated (Bernardeau 1992, Paper I) and the merits of this 
approach completely overwhelm this caveat given that the 
problem is not deterministically solved even in the simple 
Zel'dovich approximation. 

In this paper we rigorously compute the mean 9{x) 
given 8(x), i.e., (9)\s, up to third order in perturbation the- 
ory, assuming Gaussian initial conditions. As we shall see 
later, the derivation is greatly simplified when following Pa- 
per I in which the mean of 5(x) given 9{x) was calculated. 
In Paper I we explicitly calculated the coefficients, a;, ap- 
pearing in the expansion 



A(0) = (5) \ e = ai0 + a 2 {9 2 - a 2 ) + a 3 6 



(4) 



where a 2 is the variance of the velocity divergence field. The 
numerical values of dj are found to be in good agreement 
with the results of N-body simulations by Chodorowski et 
al. (1998) and by Ganon et al. (1998). The relation (§) allows 
one to reconstruct a density field from the corresponding 
velocity field. 

But why can we not invert the relation (^) to obtain 
the O(S) = {9)\s to the relevant order in perturbation the- 
ory? The reason is simply that the scatter of S around the 
mean value given by (^) introduces, in general, a bias in the 
estimate of 0(S) obtained by straightforward inversion. 



* Note a slight difference from the commonly used definition, e.| 
Bernardeau (1994) 



The paper is organized as follows: in Section ^ we de- 
scribe how the mean of the velocity divergence can be com- 
puted from the density using the results of Paper I. In Sec- 
tion ^ we calculate the numerical values of the coefficients 
entering the formula for O(S). Section ^ is devoted to the 
computation of the scatter in the 9-5 relation. Summary is 
given in Section |E| 



2 THE FORMALISM 

We now outline the derivation of the relation between the 
divergence of the velocity field and the density contrast. Let 
us express the density contrast as a sum of terms Si, each 
corresponding to the i th order in perturbation theory, 



8 = Si + S 2 + 6 3 + ■ ■ ■ , 

and, similarly, for the velocity divergence 



(5) 



(6) 



In general, the i th order solution is of the order of (Si) 1 
(Fry 1984, Goroff et al. 1986). We assume here that 9 and 
5 are well approximated by truncating the above expansion 
at the third order. The linear theory solution mentioned in 
Section 1 is simply the perturbative expansion truncated at 
the lowest, i.e. first order term, Si and 9\. (For explicit forms 
of higher-order solutions and other details see Paper I.) The 
coefficients, at, appearing in the expansion of the function A, 
equation (m, are combinations of the joint moments of the 
density and velocity divergence fields. According to Paper I 



a i 



a 2 = 



a.3 = 



1 + 

S3S 



£2 + 
- S30 



(S3S — S30) S 3 g 



£4 

2 



6 

£4 — (S3S — S 3 $) S3B 



(7) 

(8) 
(9) 



In equation ([!]), a 2 is the leading, linear, contribution 
to the variance of the density field, a 2 = (Si). The quantity 
S35 denotes the skewness of the density field, 



S38 = 



(S 3 ) _3(S?S 2 > 



(S 2 ) 2 



+ 0(a 2 ), 



(10) 



and 5*39 is the skewness of the velocity divergence field de- 
fined in an analogous way. The quantities £2 and £4 are 
given by 



£2 = 
and 

£4 = 



{5 2 9 2 ) c -{9 2 2 ) c + (8 1 53) c -(9 1 93) c 



3(8 2 S 2 9 2 ) C 



3(6>?0l) c +(SiS 3 } c - 



9'i9 3 ) c 



(11) 



(12) 



In the expressions above, the symbol (-) c stands for the con- 
nected (reduced) part of the moments. 

In the derivation of the coefficients at, no assumption 
was made about the particular forms of 5j and 9j, except 
the linear theory result Si =61. Consequently, the inverse 
relation for the mean of 9 given S, can immediately be writ- 
ten by exchanging the symbols S and 9 in equations (Q) 
and (Q)-(|l^). Thus if we express the function O = (9)\s in 
the following form 
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Q(S) = riS + r 2 (8 2 -a 2 e ) + r 3 8 3 , 



(13) 



where a 2 = (8 2 ) is the variance of the density field, then the 
coefficients n can be written as 



n 



r-2 



1 + 



E 2 + 



(S 3 



S35) S 3 s 



2 



T3 = 



S 3 e — S 3 s 
6 ' 

E 4 — (S 3 e — Sjjj) $33 



where 



E 2 



and 



E 4 = 



<0 2 5 2 ) c - <<5i)c + (6163)0 - (Si5 3 )c 



3{8 2 8 2 S 2 ) C - 3<<5 2 <5|) c + (8l& 3 ) c - (5\5z) c 



(14) 

(15) 
(16) 

(17) 

(18) 



Equations (|l^)-(|l^) are, strictly speaking, valid only 
up to third order in as and for typical values of 8, i.e., 
away from rare peaks. The spatial average of 6 is given by 
f(8)\ s p{S)dS = jQ(8)p(8)dS. The mean of the first two 
terms in equation (^) is exactly zero. The mean of the third 
term is 7-3 (8 3 ) = r 3 S 3 s&x ~ 0(aj) (see eq. Jul). Hence, our 



third-order formula ([13J) fulfills the requirement that the av- 
erage of the velocity divergence vanishes up to the terms 
cubic in as, as expected. 

By comparing the coefficients r» with the expressions 
([?])- ([[2]) for the coefficients Oj, we see that the only new 
quantities to be computed are (8 2 ) c and (8 2 S 2 ) C , appear- 
ing in E 2 and E 4 respectively. Therefore, one can express r, 
as combinations of at and some residuals Di. We uniquely 
define the residuals Di in the following way. Let us approx- 
imate the function by the perturbative inversion of the 
function A, equation (Q). The approximations of the coeffi- 
cients Ti obtained in this manner will be denoted by m. The 
exact values of the coefficients are then given by a sum of 
the approximate values and the residuals Di, 

n^m+Di. (19) 

Straightforward calculation of the coefficients rii yields 

2-a x -2a\a 2 , (20) 

-aa, (21) 
n 3 = -a z + 2at. (22) 



m 
Hi 



We will now transform the coefficients rj into forms sim- 
ilar to the coefficients rii. The case of r 2 is trivial: 

S39 — S-js S38 — S 3 e /n „v 
r 2 = g = g = -02 , (23) 

so comparing with equation (^l|) we see that D 2 = 0. The 
calculation of n and r3 is more lengthy but straightforward. 
The transformation of the coefficient r$ is simpler when we 
rewrite E 4 in a form consisting of explicitly asymmetric and 
symmetric parts, 

Sis-Su 3<6»?(<5 2 -6> 2 ) 2 } c 



v' 

^ 4 4 2a* 

Here, is the kurtosis of the density field, 

6(^ 2 2 ) c + 4(5^ 3 ) c 



Sis = 



+ <D(a z 



(24) 



(25) 



and Sig is the kurtosis of the velocity divergence field defined 
in an analogous way. The final result is 



ri = 2 — ai — 2a 2 a 2 + D\a 2 . 

r 2 = -a 2 , 

»~3 = -«3 + 2a 2 + D 3 , 



where 
Di = 
and 
D, = 



3{6i(5 2 



h) 2 ) c 



((5 2 -e 2 ) 2 ) c 5(S 3S -S 3e ) 2 



2a a 



18 



(26) 
(27) 
(28) 



(29) 



(S3 



S 3B ) 2 {8 2 (8 2 - 8 2 ) 2 ) c 



9 



2a 6 



(30) 



Note that we have introduced D 1 = D\ja 2 since this 
rescaled parameter, like D 3 , does not depend on a. 

The expressions for A and B up to second order in per- 
turbation theory can be derived by neglecting terms of order 
a A . The second order expression for the mean 8 given 8 is 
therefore 



(8)\s = 8 - a 2 {8 2 



2\ 



(31) 



(cf. Bernardeau 1992). This expression is identical to that 
obtained by direct inversion of the second order expression 
for (8)\e. Only when third order corrections are included, 
does the scatter around the mean introduce a bias in the 
relation when derived by direct inversion, ft is interesting to 
note that although Di and D 3 introduce 0(a 3 ) corrections, 
they themselves are constructed exclusively from the first 
and the second order solutions, 81 =61, S 2 and 8 2 . 



3 NUMERICAL CALCULATIONS OF n 

All of the numerical results presented in this section are per- 
formed for fields smoothed with a Gaussian filtering window 
of width R. All the terms appearing in the expressions for n 
have exactly the same mathematical structure as the terms 
appearing in the expressions for ai. Therefore, the same ar- 
guments that we used in Paper I to argue that the parame- 
ters ai are almost ^-independent apply as well to the case of 
ri. For example, r 2 = —a 2 , and in Paper 1 we explicitly de- 
rived the f2-dependence of a 2 and showed it to be extremely 
weak. Thus, the following calculations of Ti axe performed 
for a flat Universe, but we expect ri to be robust to the value 
of Q. 

The calculation of the terms {8 2 )/a 4 and (8 2 8 2 )/a e 
which contribute, respectively, to E 2 and E 4 can be found 
in Lokas et al. (1995, 1996). The rest of the terms appearing 
in the expressions for n were calculated in Paper I. 



3.1 Power law spectra 

First we present results for the scale-free spectra of the form 

P(Jfc) = Ck n , -3 < n < 1. (32) 

The coefficients r 2 and r 3 are then independent of the nor- 
malization of the spectrum and the smoothing scale. The re- 
sults are straightforward to obtain from equations (|l5|)-(|lrj) 
and they are given in Table |l| for various spectral indices. 

By comparison of r 2 and r 3 given by equations (|27|)-(|2^) 
with the approximate values n 2 and n 3 given by equations 



© 0000 RAS, MNRAS 000, 000-000 



4 Chodorowski, Lokas, Polio and Nusser 



Table 1. The coefficients n, r 2 , rz and 713 as functions of 
the spectral index n for scale-free power spectra and Gaussian 
smoothing 



n 




n 


T'2 


rz 


ri3 


-3.0 




1 + 0.3a 2 


-0.190 


0.0826 


0.0826 


-2.5 


1 


+ 0.202<r 2 


-0.192 


0.0822 


0.0822 


-2.0 


1 


+ 0.077cr 2 


-0.196 


0.0818 


0.0821 


-1.5 


1 


- 0.296cr 2 


-0.203 


0.0812 


0.0822 


-1.0 






-0.213 


0.0806 


0.0835 


-0.5 






-0.227 


0.0797 


0.0865 









-0.246 


0.0783 


0.0928 


0.5 






-0.270 


0.0756 


0.1051 


1.0 






-0.301 


0.0707 


0.1283 



(|l|)-(|2|) we find that for n > -3, the residual D 2 = but 
D3 7^ 0. Only for n = —3 we find that rz exactly equals to 713 
and D3 = 0. We present 713 in the last column of Table ^ We 
see that the approximate values 713 diverge from the exact 
ones, rz, more significantly for higher spectral indices. 

In order to obtain ri we use equation (|l4|) and get for 
the unsmoothed fields 

^ = -^-h{n)*-0A (33) 

where h(n) is the part weakly dependent on n which con- 
tributes roughly 10% to the value of T,' 2 . Using this result 
we estimate the value of the coefficient n in this case to be 

n^l + 0.3cr 2 . (34) 

The value of the residual D\ = D\ / a 2 given by equation ( |29| ) 
is for the unsmoothed fields independent of n and equal to 
-32/2205 = -0.0145. 

When smoothing is applied we are restricted to the 
spectral indices n < —1 as in the case of ai (see Paper I). 
The case of n — —3 with smoothing corresponds to the same 
case with no smoothing (because for such a spectrum the 
dominant contribution comes from the small wave-numbers 
at which the window function equals to unity), hence n is 
given by equation (|34|). For n — —2 the calculations can be 
performed analytically and we find 

29 (35) 



E' 2 



196 



-0.465 



hence 



n=l + 0.0768(7 



(36) 



For half-integer values of spectral indices we calculated 
the values of ri numerically. The results are given in the 
second column of Table [lj Comparing equation (|5(]) with 
equation @ we obtain the values of £>i = D l /a 2 . We find 
that in this range of spectral indices the coefficient m is a 
good approximation of n that is the correction introduced 
by Di is very small. Therefore, in Table |l| we do not present 
the values of m. 



3.2 The standard CDM spectrum 

As an example of a scale-dependent power spectrum we 
consider the standard CDM (with coefficients given by Efs- 
tathiou, Bond & White 1992) normalized so that the linear 
rms fluctuation in spheres of radius R = 8/1 -1 Mpc is equal 



-0.95 
1.2 



1.0 
0.8 
0.6 
0.4 
0.2 
0.0 
-0.2 
-0.4 



-0.52 



-0.26 



-0.07 



LUM 


ri 













10 15 

R (h _1 Mpc) 



20 



Figure 1. The coefficients n, r 2 and r$ for the standard CDM 
spectrum normalized to (top hat) u% = 1 in the weakly nonlinear 
range of Gaussian smoothing scales. The values of the effective 
index corresponding to each of the smoothing scales are also plot- 
ted. 



to unity. Given w> the coefficient r 2 is obtained immediately 
from equation (|27|). 

In calculating rz we use the effective index 



n c ff = 



R da 2 (R) 
a 2 dR 



(37) 



which measures the slope of the power spectrum at smooth- 
ing scale R. Then the values of rz at a given scale are found 
by interpolating the values known for scale-free spectra (Ta- 
ble |l]) at the effective index corresponding to that scale. 

The values of T,' 2 that contribute to ri cannot however 
be calculated using the effective index method. We therefore 
calculate the correction to unity in n by numerical integra- 
tion of the relevant formulae (see Paper I) . As in the case 
of ai we benefit from the fact that the slope of the CDM 
spectrum approaches —3 at large wave-numbers and thus 
the weakly nonlinear correction to the value of n converges 
even for scales where the effective spectral index is close to 
+1. 

The results for the CDM spectrum in the weakly non- 
linear range of scales are presented in Figure |l| In the values 
of n we incorporated the linear variances a 2 determined by 
the normalization to as = 1. The Figure also shows the 
values of the effective index corresponding to each of the 
smoothing scales. 



4 SCATTER IN THE 8-6 RELATION 

In Paper I we have computed mean & given 8, {5)\e- Analo- 
gous calculations lead to a formula for the conditional vari- 
ance, a 2 \e = ((6 — (5)|e) 2 )|e; we present them in brief in 
Appendix A. The result is 

— b ag + b 2 cr%6 2 + O(o-f) , 



a 

where 



60 = 



((8 2 -e 2 ) 2 ) (6 2 (5 2 



h) 2 ) , (ASz) 2 



+ 



(38) 



(39) 



2<r 6 18 
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Table 2. The coefficients f>o and b 2 as functions of the spectral 
index n for scale-free power spectra and Gaussian smoothing 

n b b 2 



-3.0 
-2.5 
-2.0 
-1.5 
-1.0 
-0.5 
0.0 
0.5 
1.0 

and 



32 

2205 



! 0.0145 

0.0155 

0.0193 

0.0288 

0.0517 

0.116 

0.378 



3.61 
2.66 
1.02 
2.86 
6.77 
1.45 
2.95 
5.76 



(A5 3 



!) 



(40) 



2cr6 

Here, A S3 = Sas - S 3 e- 

In the present paper we have computed mean 8 given 
8, (6)\s- The conditional variance in this case, a 2 \s = ((8 — 
(0)\s) 2 )\s, may be similarly obtained from equations (Bq)- 
(kOJ) by exchanging the symbols 8 with 8. Unlike the coeffi- 
cients of the mean trend, equations (f^)-(p^), the coefficients 
60 and 62 are invariants with respect to this exchange. There- 
fore, 



a 2 \ s = bo(T$ + b 2 a 2 5 2 + C(cr|) . 



(41) 

and @ respectively. 



with 60 and 62 given by equations (|3 

The coefficients bo and 62 may be written in terms of 
the residuals D\ and D3, equations (E9f) and rt3C| 



6 = -£>i 
and 

62 = -D, 



2D 3 



(42) 



(43) 



Thus the presence of the residuals in equations fl2q)- 
( p8^ for the coefficients is a direct consequence of the 
scatter in the relation between 5 and 8. Were the relation 
entirely deterministic, it would be describable by one func- 



tion, so then O = A 



and there would be no resid- 



uals. However, since in weakly nonlinear regime the values 
of pairs (5,6) form an elongated set of some scatter, it is not 
surprising that averaging along different coordinates gives 
different curves. 

As typically 5 ~ 6 ~ a, we have a 2 8 2 ~ a 2 8 2 ~ a 4 
and both terms entering the formulas (^i|) or ( [il"| ) for the 
conditional variances are of the order of a 4 . The rms value 
of the scatter in the 5-8 relation - the square root of the 
conditional variance - is therefore of the order of a 2 . A good 
measure of the elongation of a two-dimensional set of points 
(5,8) is the ratio of its scatter to a typical value of 5 or 8. 
In our case this ratio is ~ a. When a — > 0, the ratio thus 
tends to zero, as expected: in the limit of linear theory the 
relation between 5 and 8 is deterministic. 

We present the values of bo and 62 as functions of the 
spectral index n in Table ^| We see that bo and &2 are always 
positive and the conditional variance is positive-definite, as 
required. 

It is rather astonishing that the formulas for the coeffi- 
cients bo and i>2 for the conditional variance are constructed 
exclusively from second-order perturbative contributions for 
5 and 8. In section 2 we have shown that up to second order 



Q = A -1 ; the residuals Di appear only when third-order 
corrections to the mean trend are included. Still, Di's and 
bi's contain only second-order terms. Why? The explana- 
tion of this apparent paradox is provided by Chodorowski 
(1997). In this work, for unsmoothed fields, the weakly non- 
linear density has been shown to be a local function of the 
two velocity scalars: the expansion (velocity divergence) and 
the shear, 



5(x) = 8(x) + 



21 



2 (x)-^ 2 (x) 



0(a 3 



(44) 



(cf. Mancinelli & Yahil 1995 and Catelan et al. 1995). Here, 
the shear scalar E is 



E = (EijEij) ^ , 

where 

_ 1 1 

and Vij are velocity derivatives. Thus, the scatter in the 8— 
8 relation comes from the shear term in equation (fl4|). The 
mean trend is 



(45) 



(46) 



(Si\e = 0+21 



+ 0(a 3 



(47) 



It reproduces up to second order equation (||) since at the 
lowest order 8 and E are statistically independent, hence 
(E 2 )|e = (E 2 ) = (2/3)cr 2 (see Chodorowski 1997 for details). 
The conditional variance is 



4 

49 



49 V 



(E 2 )) : 



0(a :i )Y 
+ 0(a 5 ) 



(48) 



The linear term in equations ( f44| ) and (|47j), 8, has can- 
celled out in the above formula for the conditional variance. 
Therefore, third-order contributions to equation (^ij) result 
in contributions to the conditional variance that are already 
of the order a 5 . This is in contrast to the formula for weakly 
nonlinear growth of ordinary variance, 



(S 2 )-(5 2 ) = (5 2 )+ 2(5 1 S 3 ) 



(49) 



which does involve a third-order contribution due to the 
presence of the linear term. 

The above consideration is strictly valid only for un- 
smoothed fields, but the linear terms in the conditional vari- 
ance will also cancel out for smoothed fields and the lowest 
order contribution to it will be constructed exclusively from 
second-order terms. 

The ensemble average on the RHS of equation (^8|) can 
be performed: the result is 



32 4.^/5 

2205 V 



(50) 



The case of unsmoothed fields corresponds to the case of 
smoothed fields with the spectral index n = — 3. From Ta- 
ble |2| we see that for n = —3 the parameters 60 = 32/2205, 
02 = and equation ( p8[ ) coincides with the above formula. 
Thus, equation (^) is an extension of the result of Chodor- 
owski (1997) for the case of smoothed fields. 

In this work we are mostly interested in the conditional 
variance {(8 — {8)\s) 2 )\s- In Table [| we see that b 2 <C bo al- 
ways. For mildly nonlinear fields (as ^ 1) this implies that 
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in equation (^l|) for <t 2 |<s the term quadratic in 5 is negli- 
gible with respect to the constant one. [The coefficient 62 
is exactly zero in the case of a top-hat smoothing (Chodor- 
owski et al. 1998)]. Furthermore, the actual 5-dependence 
of the conditional variance must originate from the terms 
that are of order higher than a\. In Appendix A we show 
that the next-order term is tx aj5 which is of the order of 
erf; other terms are already O(af). This is indeed observed 
in N-body simulations, where the estimated variance can be 
well described by the formula (Chodorowski et al. 1998) 

2\(NB) ,(NB) 4 . ,(NB) 4; 

a \s = b o °~s + b i °~s5, (51) 

with b^" 8 ^ ~ O.Sbg NB ' 1 . In principle, the calculation 
of the coefficient 61 involves fourth-order perturbation the- 
ory. If density is a local function of the velocity diver- 
gence and the shear up to third order, similar arguments 
to these presented after equation (^) may be used to argue 
that the formula for 61 should be composed from at most 
third-order perturbative contributions. Even if this is true, 
however, we find that the moments a~ 6 {S2(S2 — #2) 2 )} and 
a~ 6 (9i(S2 — #2X^3 — #3)} will enter into this formula, of extra 
complexity in comparison to any moments computed here. 
Given all that, in this paper we will not try to predict the 
value of bi. 



the Tully-Fisher relation). As an example, let us consider 
the /R^IS-potent comparison (Sigad et al. 1998), employ- 
ing Gaussian smoothing length of 12 h~ 1 Mpc. At this scale, 
ag is approximately 0.3, hence the ratio of the rms value 
of the scatter around our relation to the rms value of the 
divergence lies in the range 0.04-0.18. In contrast, the rms 
value of the observational scatter is approximately 0.21 (es- 
timated from the mock catalogs in the 'standard volume'), 
thus the observed ratio is as high as 0.7. N-body simula- 
tions (Nusser et al. 1991, Mancinelli et al. 1994, Chodor- 
owski et al. 1998; cf. also Sigad et al. 1998) also demonstrate 
that the scatter around relations of this type is relatively 
small. The relation and the formalism presented here can 
be valuable for quantifying the redshift space anisotropies 
in the non-linear regime. This is particularly important for 
probing small scales in the planned SDSS and 2dF surveys. 
We intend to use the relation for estimating the amount of 
anisotropy of structure in redshift space. 



ACKNOWLEDGMENTS 

This research has been supported in part by the Polish State 
Committee for Scientific Research grants No. 2.P03D.008.13 
and 2.P03D.004.13. 



5 SUMMARY 

We have derived a local relation between the divergence 
of the weakly non-linear peculiar velocity field and the 
corresponding density contrast field. Specifically, we have 
computed the mean value of the velocity divergence given 
density up to third order in perturbation theory, assum- 
ing Gaussian initial conditions. Our perturbative calculation 
yields a third-order 'expansion' of the mean value of the di- 
vergence in terms of the density contrast. The coefficients of 
this expansion have been explicitly calculated for scale-free 
and standard CDM power spectra, for Gaussian smooth- 
ing of the fields in question. In the case of CDM, the co- 
efficients depend weakly on the smoothing scale. Moreover, 
they are expected to depend very weakly on Q. It is inter- 
esting that the value of the linear coefficient in the relation 
differs from unity - the value predicted in linear theory. The 
corrective term to this value is proportional to the variance 
of the density field. The same form of the corrective term is 
found in the relation between the mean velocity divergence 
and the density resulting from the Zel'dovich approximation 
(Chodorowski 1998). 

In order to assess the tightness of the relation, we have 
computed the scatter in the velocity divergence around its 
mean value given density. We find that, at least up to third 
order in perturbation theory, the scatter is almost indepen- 
dent of the density contrast. The rms value of the scatter rel- 
ative to the rms value, ag, of the divergence is approximately 
1/2 

b ag where bo — 0.015-0.37 for scale free power spectra 
with n ranging from —3 to 0. Over scales larger than about 
5ft _1 Mpc, this scatter is small relative to that introduced 
by observational errors in analyses of redshift surveys and 
catalogs of peculiar velocities, such as the density-density 
or velocity- velocity comparisons. The main source of these 
errors are the sparse and non-uniform sampling of the ve- 
locity data and the scatter in the distance indicators (e.g. 
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APPENDIX A: CALCULATION OF THE CONDITIONAL VARIANCE 

Here we describe the calculation of the variance a 2 \g — (S 2 }\o — (S)\g. The derivation is analogous to the derivation of the 
mean trend, presented in Paper I, therefore we will sketch it in brief. If p(S, 9) is the joint probability distribution function 
(PDF) for 5 and 9, then mean 5 2 given 9 is 



(S 2 )\o = 



f 5 2 p(S,9) dS 
W) 



The quantity N = J 5 2 p{S, 8) d6 can be expressed as 



N - 



1 

2^ 



e isB -JL^^is) 



ds . 



(Al) 



(A2) 



where $ is the characteristic function of the joint PDF. It is related to the cumulant generating function, K,, by the equation 
$(it, is) = exp [K(it, is)] . (A3) 
The cumulants, K mn , from which K, is constructed, 



K. = 



E 



r(W) m (i*r 



(m,ra)^(0,0) 

are given by the connected part of the joint moments 

Kmn = {8 9 }c • 



Using equations (A3) and (A.4) we obtain 

d 2 



d(it) 



Mit,is) 



EK2n ,. % n . \ * KlnKlp ,. xn+p 
— ( is ) + 2^ -r^i-W 

n — n,p=0 



nip] 



exp 



^r (ls) 

n=l 



By defining z — k^ 2 s, /i = 6/crs = S/k^q 2 , v = 9/ag = O/k^ 2 , and the standard cumulants A mn = K m n/(«: m/ ' 2 '" n/ ' 2 
that 



(A4) 



(A5) 



(A6) 



20 K 02 J, we find 



Af = 



1 K 2 Q 

2ty 1/2 



dze z { 



n\ ^ nip! 

.ro— n,p— 1 



exp 



n=3 



(A7) 



In weakly nonlinear regime the standard cumulants obey the following scaling hierarchy (Fry 1984, Bernardeau 1992) 



m + n — 2 



(A8) 



where a is the linear variance of 5 or, equivalently, of 9 (recall that at linear order 5 = 9). The series in equation (A7) are 
thus power series in a small parameter a and we can truncate them at some order p neglecting contributions which are of 
order > a p . The leading-order formula for the conditional variance is obtained by keeping the terms up to the order of a 2 . 
Integrating the resulting expression yields 

N = 1 



«20 -i„ 2 

e 2 x 



27r k]1 2 



H {v) + \\iHi(u) + Aai-ffiM + ( 



Al2 + 



6 



H 3 {u) + ^-H B (v)+ 
6 



-H 2 (v) + 



Al3 

3 



24 



Ao4 . A12A03 . A03 



24 



+ 



72 



+ 



H 6 (u) + ^Hs(v) 



Hi{v) + 



(A9) 



where H n are the n-th order Hermite polynomials. 

We now turn to calculating the PDF p(9) in equation (|Al[ ). Computing it in a similar way as M we rederive the so-called 
Edgeworth expansion for the variable v (Longuet-Higgins 1963, Bernardeau & Kofman 1995, Juszkiewicz et al. 1995), 



From equations 



-1/^/2 



1 + gAosfliW + — \ qa Ha{v) +^\i 3 H a {v) 



and (A10), after some algebra we obtain 
where 



(A10) 



(All) 
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So(u) 



1 - Afx + Ani/ 2 . 



Si(u) = (A21 - 3Ai2 + 2A 03 )^ + (A12 - A 3)^ 3 , 
and 



(A12) 
(A13) 



S 2 {u) = 



A22 — 2A13 + A04 . 2A21A03 + 3Ai2 — IOA12A03 + 5A§3 



+ 
z 

4Ai3 + 3Ao4 — A21A03 



+ 



3Af 2 + I2A12A03 — 8A03 



A13 — A04 A?2 — 6A12A03 + 5Aq 3 



(A14) 



Note that Sk{v) are of the order of o h . To obtain (/i 2 )|i/ — (a')Ik we use the expression for given explicitly in Paper I. 

The result is 



(^ 2 )\ v -{^)\i = 1 - An 



A22 — 2A13 + A04 , A21A03 + A 12 — 4A12A03 + 2A 



+ 



(A 2 i - 2A12 + A 03 )^ + 
A22 — 2A13 + A04 



2A13 + Aq4 A21A03 + 2A12 — 6A12A03 + 3Aq3 



(A15) 



Were the variables identical (to all orders), \i = v, then the cumulants would be An = 1, A mn = A( m+n ) , and the conditional 
variance would be zero, as expected. Returning to the 'physical' variables S and 8, using the leading-order expressions for the 
standard cumulants, like 

{(s 2 - e 2 f) 



1 



and recalling that ag = a + 0(a 3 ) we finally obtain 

<«5 2 )| fl - <5)| e 2 = boat + b 2 a 2 g 6 2 + 0(a 5 g) , 
where 

, {(s 2 -e 2 f) {el{8 2 -e 2 f)) (as 3 ) 2 

°0 = ; 7T- 7, 1- 



and 



2a 6 



(el{8 2 - e 2 ) 2 )) (as 3 ) 2 



18 



2(7 6 



9 



(A16) 
(A17) 
(A18) 

(A19) 



In equation (A15) higher-order corrections to the values of cumulants yield corrections to formula (4.17) which are already of 
the order of ag. The only exception is the term linear in v = 6 /erg which yield s a c ontribution of the order of ag. Therefore, 
for mildly nonlinear fields, the first departure from the leading-order formula ( A17 ) for the conditional variance a 2 \g should 
have the form oc agO. Analogously, the expression for the conditional variance in the reverse case, a 2 \$ = (6 2 }\s — should 
contain a higher-order term of the form oc ajS. 



This paper has been produced using the Royal Astronomical Society/Blackwell Science I^TffjX style file. 
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